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Abstract 



Quasinormal mode spectra for gravitational perturbations of black holes in 
four dimensional de Sitter and anti-de Sitter space are investigated. The de 
Sitter case lends itself to approximation by a Poshl- Teller potential. The 
anti-de Sitter case is relevant to the ADS-CFT correspondence in superstring 
theory. The ADS-CFT correspondence suggests a prefered set of boundary 
conditions. 
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I. INTRODUCTION 



The spectrum of gravitational waves emitted by a black hole which has been disturbed is 
dominated, over intermediate timescales, by characteristic frequencies and decay timescales. 
These can be accurately described in terms of the discrete set of quasinormal modes. The 
frequencies of these modes can be determined numerically these days to high accuracy up 
to large order 

A new application of the quasinormal mode spectrum has arisen recently from superstring 
theory. There is a suggestion, known as the ADS-CFT correspondence, that string theory in 
anti-de Sitter space (usually with extra internal dimensions) is equivalent to conformal field 
theory in one less dimension A black hole, which has a characteristic temperature fixed 
by the Hawking effect, should correspond to CFT at finite temperature and the quasinormal 
modes would describe non-equilibrium effects @. 

The temperature of the black hole is fixed by the surface gravity of the event horizon 
Ki. The anti-de Sitter space also has a characteristic length scale a which is conventionally 
refered to as the anti-de Sitter radius. In the non-rotating case, there are two black hole 
solutions for a given value of the surface gravity |0, a larger one with positive and smaller 
one with negative specific heat. Therefore the metric is determined by a, ki and the choice 
of the larger or the smaller black hole. 

Previous work on the quasinormal mode spectrum in anti-de Sitter space has concen- 
trated on the wave equation for scalar fields. Chan et al. studied the scalar field problem 
prior to the interest in the ADS-CFT correspondence. Horowitz et al. P,p!Q| examined the 



scalar spectrum in four and five dimensions. They found simple scaling behaviour of the 
spectrum in both limits ri ^ a and ri -C a, the former explained by an approximate 
symmetry of the metric in that limit. 

Gravitational perturbations of the BTZ black hole in three dimensions and the anti-de 
Sitter black hole in four dimensions have been studied by Cardoso et al. [0,0- They make 
use of the fact that the gravitational perturbation equation is separable in four dimensions 
T3|-15 . This fact was also used in a preliminary report on our results [1^. The main 



difference between our work and the work of Cardoso et al. is in the choice of boundary 
conditions on the axial and polar metric perturbations. 

There is a symmetry which relates the axial and polar metric perturbations, whose 
existence can be explained by the fact that the covariant form of perturbation analysis 



based on the Weyl tensor is governed by a single equation |]18[. We argue that the thermal 
fluctuations in the energy of the conformal field theory are determined by a single set of 
modes for which the axial perturbations have dirichlet boundary conditions and the duality 
between axial and polar perturbations is preserved. We calculate the quasinormal mode 
frequencies numerically. For fluctuations in other components of the stress energy tensor 
there is a one parameter family of boundary conditions and the modes depend weakly on 
this parameter. 

We have also investigated the quasinormal mode spectrum for black holes in de Sitter 
space. An interesting feature of this problem is that the asymptotic behaviour of the modes 
is under better control than the corresponding problem in flat space |l| . The first few modes 
where evaluated previously in ref. [jl3[. Using methods based on continued fractions now 
gives us many thousands of modes. A Poshl- Teller approximation for the potential works 
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extremely well in this example . 



II. DE SITTER SPACE 

The spherically symmetric de Sitter spacetime which contains a black hole has two hori- 
zons, one at the event horizon of the black hole and a cosmological horizon further out. In 
standard coordinates, the metric for a black hole of mass M and for a cosmological constant 
A can be written 

ds"" = -^df + ^dr^ + r\de^ + sin2 6 #2) (i) 

where 

A = - 2Mr - ^Ar^ (2) 

The event horizon lies at the smallest root ri of A, where the surface gravity Ki = r^^A'/2. 
We can invert these relations to give 

M = |riKri + l) (3) 
A = rf{l - 2Kin), (4) 

showing, incidentally, that Kiri < 1/2 in de Sitter space. 

The gravitational wave equations were separated originally for charged black holes [jl3 



then for uncharged and finally for rotating black holes [T^. The metric perturbations can 



be classified by their axial or polar symmetry under spatial inversions. The transformation 
theory relating the two types is described in the next section. 

The axial metric perturbations with time dependence e*'^* and angular mode / are subject 
to a one dimensional scattering equation. 



d^ij 



dr 

where 



5 + p-V(r)}^. = (5) 



Vir)^^(Hl^l)-'-^) (6) 



^4 



The tortoise coordinate r*, given as usual by dr* = dr / A, covers the range — oo < r* < oo 
as r ranges between the two horizons. The potential decays to zero exponentially fast in 
the r* coordinate at both horizons. This is an important property which is not present for 
the black hole in an asymptotically flat spacetime and guarantees that the reflection and 
transmission coefficients are meromorphic functions of the frequency . 



By definition, the quasinormal modes are ingoing waves at the event horizon and outgoing 
waves at the cosmological horizon. 



^ . * (7) 

' "^^'^ as r* ^ -oo ^ ' 
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The method which we use to determine the complex frequencies follows the continued frac- 



tion approach devised by Leaver The idea is to first scale out the divergent behavior 
at the cosmological horizon and examine the resulting Frobenius series. If we begin with a 
purely divergent solution at the event horizon, the values of uj can be chosen to select the 
convergent rescaled solution at the cosmological horizon. 

The analysis is slightly easier if we use the independent variable x = 1/r. In de Sit- 
ter space, the radial equation has three regular singular points xi (the event horizon), 
X2 (the cosmological horizon) and x^. The indices of the equation at these points are 
±Pj = ±ia;/(2Kj), where 

K^ = 2M{xi-X2){xi-x^) (8) 

K2 = 2M{X2 - X'i){x2 - Xi) (9) 
/€3 = 2M(X3-Xi)(x3-X2) (10) 

Note that pi + P2 + Ps = and 

e-'^^' = (x - xiyp' (x - X2Y'"' (x - x^)-P^ (11) 

We define, 

^ = e-^^"*M (12) 

For a quasinormal mode, u{x) must be divergent at xi and finite at X2- The Frobenius series 
is therefore 

u{x) = (x - x,fP^ a„ (l^l^y (13) 

„=i \X2-X1J 

where Qq = 1 and the coefficients satisfy a recurrence relation of the form 

anttn+l + Pnttn + Indn-l = 0. (14) 

The existence of this recurrence relation with only three terms is very important, because 
the convergence of the Frobenius series at x = X2 is now dependent on a continued fraction 
equation, 

which can be solved very quickly numerically. 

Figures |l| and ^ show the real and imaginary parts of the frequency of the I = 2 and I = 3 
modes as ri^i is varied. The quasinormal modes in flat space are recovered at ri^i = 0.5. 

The simple shape of the potential suggests that a Poschl- Teller approximation might be 



effective, as it is for the Schwarzchild black hole |2^. The approximate potential 



VpT = Vosech'{r*/b) (16) 

contains two free parameters which are used to fit the height and the second derivative of 
the potential at the maximum. The location of the maximum has to be found numerically. 
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The quasinormal modes of the Poschl- Teller potential can be evaluated analytically, 



a;„ = ^ |(n + l/2)z + ^fenW/i} . (17) 

for n = 1, 2, . . .. These values are compared to the numerical values for Z = 2 in figure |l| and 
/ = 3 in figure 0. The approximation works best for the imaginary parts of the frequencies, 
but the differences between the numerical values and the Poschl- Teller values decrease as / 
is increased and soon become smaller than the resolution of such graphs. 

For large / it is possible to obtain the Poschl- Teller fit to the potential analytically. The 
maximum of the potential in this case lies close to r = 3M. From the second derivative at 
r = 3M, we get 

= /ci(l + Kin)-\l + |/€iri)i/2 (18) 

and 6^Vo = /(/ + 1). Equations and (^ have been used to express the result in terms of 
the surface gravity. The analytic expression for the frequencies is therefore 

= K^-^^i^^^ Un + 1/2)^ + (/(/ + 1) - l/4)Vn . (19) 

Figures ^and Q show the numerical results and the approximation for small /. The difference 
is largest in the fiat space limit at the right hand side of the figures. As / is increased, the 
difference beween the approximation and the numerical results soon becomes insignificant. 



III. ANTI-DE SITTER SPACE 

The metric for a black hole in anti-de Sitter space takes the same form as for de Sitter 
space but with A negative. The cosmological horizon is no longer present. We introduce an 
anti-de Sitter radius a by the definition = —3/ A. For fixed values of Ki and a there are 
two black hole solutions with different radii, 

1/2 



Kin = lii^ia)^ ± i [{Kiaf - 3{Kiaf) (20) 



The larger solution has positive and the smaller one negative specific heat. 

The axial metric perturbations are subject to the same radial equation as for de Sitter 
space, however in anti-de Sitter space the tortoise coordinate lies in the range — oo < r* < 0. 
The boundary conditions at r* = can be dirichlet, neumann or robin depending respectively 
on whether the field, its derivative or a combination of both vanishes. Previous work has 
tended to opt for dirichlet boundary conditions [Q, based largely on ADS invariance scalar 



field theory p^ , |25[| . However, dirichlet boundary conditions break a duality in the equations 



between axial and polar metric perturbations |T^. We will introduce a set of self-dual 
boundary conditions based on the ADS-CFT correspondence. 

In the Newman-Penrose formalism, the gravitational perturbations can be described by 
two Weyl scalars \l/o and \l/4. In flat space, \I'o would represent gravitational waves which are 
incoming at infinity and \I^4 would represent outgoing gravitational waves. The equations 
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for the Weyl scalars separate to give a generahsation of the Teukolski equations |T^. We 
put, for example, 



3 

*o = e^"*^r(r)5,(e) (21) 
where Si{6) can be found in Chandrasekhar |T^. The radial equation for Y{r) can be related 



to the radial equation for the metric perturbations if we set 

Y = V^^^ + {W^ + 2iiu)A+^p^ (22) 

where 

= A + (23) 

There are two equivalent formulations corresponding to polar and axial metric perturbations, 
denoted by the plus and minus signs. The axial form, corresponding to has V~{r) = V{r) 
given previously (||). 

The axial and polar functions are related by 

C^^T = (C^ + 72M^f)ilj^ T 12MK_^^ (24) 

with C± = /i2(;x2 _^ 2) ± l2iMuj, and 

^ ^ ^Ji^TToM' ^^^^ 

where /x^ = (/_!)(/ + 2). 

The ADS-CFT correspondence conjectures that when the Anti-de Sitter space is confor- 
mally compactified, the generating functions of the conformal field theory on the boudary 
are related to the partition function of the superstring theory on the interior, or 

exp (^J h.^T^''^ ^^^^ = Z (h,,) (26) 

If we calculate the correlation fuctions of a polar quantity, such as the energy, we must 
therefore set the axial metric perturbation to vanish on the boundary. The transformation 
theory (p^ ) then implies a Robin boundary condition on the polar metric perturbation, i.e. 
the normal derivative is determined by the value of the field. The same set of quasinormal 
modes is obtained from the axial and polar equations and the duality between the two is 
preserved. 

In general, both the axial and polar metric variations are non-vanishing. This is consis- 
tent with the conformal theory viewpoint, where the six components of the metric perturba- 
tions in three dimensions are subject to three coordinate transformations and the conformal 
symmetry, leaving two independent functions. If we define the ratio 7 = , the two 

metric variations are related by (|2^), which gives a robin boundary condition 

(e + A-)r (27) 
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where 

C+ 6M 
12M ^ ^ 



e=(l-7)^ + ^ (28) 



Some special solutions lead to vanishing of the Newmann-Penrose scalars and might be 
regarded as pure gauge modes from the point of view of the ADS theory. For example, 

(29) 

with uj = —jj^i^jj? + 2)i/12M, obtained by setting ip^ = 0, is a purely axial perturbation 
with y = 0. Replacing uj by —uj gives a quasinormal mode with frequency 

Such modes appear in the numerical solutions discussed below. 

The frequencies of the quasinormal modes can be found numerically by adapting the 
method used for de Sitter space. We first scale out the behaviour of the quasinormal mode 
at the event horizon, 

^- = e'""*M, (31) 

where u will be regarded as a function of x = 1/r. The radial equation (||) gives the equation 

pu" +p'u' + 2iuj u' + {/(/ + 1) - 6Mx} M = (32) 

where p = — — + 2Mx^. We take the Frobenius series which is regular at the event 
horizon, 

u=J2^n (- — (33) 



n=0 



-Xi 



with ao = 1. The coefficients are obtained from a recurrence relation similar to ([IJ)- The 
series can then be substituted into the boundary condition m = at a; = to obtain an 
equation for the frequencies of the quasinormal modes with vanishing axial perturbations. 

The frequencies of the first few modes using dirichlet boundary conditions on the axial 
perturbations have been plotted for ani = 2.0 and ri = a in figure ^ For I = 2, one of the 
frequencies lies on the imaginary axis. (This is also found by Cardoso et al We have 

found a single frequency on the imaginary axis and can trace it for KiTi > 1.3 when / = 2, 
for KiTi > 3.4 when / = 3 and for KiTi > 38.0 when I = 10. 

The dependence of the first three modes on the metric parameters, with / = 2 and 
/ = 10, is displayed in figures |^ and |[ Both the real and imaginary parts of the frequencies 
become proportional to Ki for large ki, which agrees with the behaviour for the scalar field 



quasinormal modes |]T0|- For ri -C a, the numerical data suggests that imaginary part of 
the frequency approaches zero as in the scalar case [pUf]. 

Corresponding results for dirichlet boundary conditions on the polar perturbations (7 = 
0) are shown in figures |^ and |^. The special mode in this case is the 'pure gauge' mode (|30D 
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and the analytic solution provides a good check on the accuracy of our numerical technique. 
The frequencies for these boundary conditions are not drastically different from the dirichlet 
axial boundary conditions. The two sets of results represent two extreme cases of the one 
parameter family of boundary conditions (^) given above, and indicate a weak dependence 
on the parameter 7. 

The potential V{r) for anti-de Sitter space no longer vanishes as r ^ 00 and it is harder 
to fit the potential to a Poschl- Teller form. More importantly, the value of r* at infinity is 
zero, and the asymptotic forms of the mode functions for the Poschl- Teller potential can no 
longer be used. 

IV. CONCLUSION 

The quasinormal modes for the gravitational perturbations of a black hole in de Sitter 
space have much in common with the quasinormal modes for the black hole in asymptotically 
flat spacetime. The radial scattering problem in de Sitter space is better behaved because 
the potential falls off exponentially rather than as a power law. This explains why the 
approximation to the quasinormal mode frequencies, derived in section 2 from a Poschl- 
Teller potential, is in such good agreement with the numerical results for even moderate 
angular order /. 

The corresponding problem in anti-de Sitter space requires careful consideration of the 
boundary conditions at infinity. The proposal here has been to use the ADS-CFT corre- 
spondence and the duality symmetry between axial and polar metric perturbations. For 
correlation functions of the energy in the conformal field theory, the appropriate boundary 
conditions are that the axial metric perturbations vanish and the polar perturbations satisfy 
a robin boundary condition. 

In general, two functions have to be specified on the boundary and there is a one pa- 
rameter family of boundary conditions. The effect on the quasinormal frequencies is quite 
small. 

It would be desirable to analyse thermal widths in three dimensional conformal field 
theory to look for signs of a correspondence with the quasinormal mode frequencies found 
here. This could give an impetus for analysing the five dimensional gravitational quasinormal 
modes, where a full separation of the gravitational perturbation equations is still an unsolved 
problem. 
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FIGURES 




FIG. 1. Real (left) and imaginary (right) parts of the gravitational quasinormal mode frequen- 
cies riuj with I = 2 plotted as a function of the surface gravity ri^i, together with the Poschl- Teller 
approximation. 




FIG. 2. Real (left) and imaginary (right) parts of the gravitational quasinormal mode frequen- 
cies riU! with I = 3 plotted as a function of the surface gravity ri^i, together with the Poschl- Teller 
approximation. 
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FIG. 3. The frequencies of the first few quasinormal modes for the black hole anti-de Sitter 
metric with dirichlet axial boundary conditions and ani = 2.0 (and ri = a). Results for Z = 2, 3 
and 10 are shown. 
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FIG. 4. The dependence of the first two quasinormal mode frequencies on the metric parameters 
for I = 2 is shown for the dirichlet axial boundary conditions. The frequency with lor = is also 
plotted and designated 'special'. 
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FIG. 5. The frequencies of the first few quasinormal modes for the black hole anti-de Sitter 
metric with the dirichlet polar boundary conditions (7 = 0) and Kia = 2.0 (and ri = a). Results 
for Z = 2, 3 and 10 are shown. 




FIG. 6. The dependence of the first two quasinormal mode frequencies on the metric parameters 
for Z = 2 is shown for the dirichlet polar boundary conditions (7 = 0). The frequency with lor = 
is also plotted and designated 'special'. 
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